Fractional topological phase on spatially encoded photonic qudits 
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We discuss the appearance of fractional topological phases on cyclic evolutions of entangled qudits 
encoded on photonic degrees of freedom. We show how the spatial correlations between photons 
generated by spontaneous parametric down conversion can be used to evidence the multiple topo- 
logical phases acquired by entangled qudits and the role played by the Hilbert space dimension. A 
realistic experimental proposal is presented with numerical predictions of the expected results. 
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I. INTRODUCTION 



Geometrical phases are robust means for implementing 
unitary gates useful for quantum information processing 
[l], [|J . The phase evolution of entangled qubits was in- 
vestigated in refs.Jl, and the role of entanglement on 
the topological nature of geometric phases has been dis- 
cussed both theoretical [5-7] and experimentally Q in 
the context of two-qubit systems. More recently, it has 
been shown that the dimension of the Hilbert space plays 
a crucial role on the topological phases acquired by entan- 
gled qudits [To[ • The appearance of fractional phases is a 
remarkable prop erty of two-qudit systems, also shared by 
multiple qubits However, experimental demonstra- 
tions of these interesting features are still rare. Fractional 
phases were originally investigated in quantum Hall sys- 
tems in connection with different homotopy classes in the 
configuration space of anyons. Moreover, the topological 
origin of fractional phases in quantum Hall systems open 
interesting perspectives related to fault tolerant quantum 
computation [12j |. 

In the present work we propose an experimental setup 
to observe fractional topological phases with entangled 
qudits encoded on the transverse position of quantum 
correlated photon pairs generated by spontaneous para- 
metric down conversion (SPDC) |l3l — Fl5j | - Local unitary 
operations can be applied in this degree of freedom with 
the aid of spatial light modulators that introduce a pro- 
grammable phase profile on the wavefront of the quan- 
tum correlated photons. Then, the desired fractional 
phases can be observed through polarization controlled 
two-photon interference [H|,[lj]]- The roles played by en- 
tanglement and Hilbert space dimension are unraveled 
by the behavior of the interference fringes as the unitary 
operations are continuously changed, raising the conjec- 
ture that entanglement and dimensionality witnesses are 
involved in the quantities to be measured. These ideas 
are supported by analytical and numerical calculations of 
the expected interference fringes. The paper is organized 
as follows: in section II we introduce the fractional topo- 
logical phase concept; the experimental proposal is pre- 
sented in section III; section IV is devoted to numerical 
simulations, and in section V conclusions are outlined. 



II. FRACTIONAL TOPOLOGICAL PHASES 



Let \ip) = J2i j=i a ij\ij) be the most general two- 
qudit pure state. We shall represent this state by the 
d x d matrix a whose elements are the coefficients Qy. 
With this notation the norm of the state vector becomes 
(tjj\ijj) = Tr(crQ;) = 1 anc l the scalar product between two 
states is ((j>\tp) — Tr(/3 T a), where f3 is the dxd matrix con- 
taining the coefficients of state \<f>) in a given basis. In or- 
der to characterize a general vector in the Hilbert space, 
ote that any invertible matrix admits a polar decomposi- 
tion a = e 1 ^ Q S, where Q is a positive definite Hermitian 
matrix, and S € SU(d). Following the usual terminol- 
ogy of spontaneous parametric down conversion, we shall 
label the two qudits as signal (s) and idler (i). Under lo- 
cal unitary operations U s and V, the coefficient matrix is 
transformed according to a' — U s a V^ . It can be read- 
ily seen that this kind of transformation preserves the 
polar decomposition and can be represented separately 
in each sector of the coefficient matrix: a' = e 1 ^ Q' S' , 
where e_^' = ^det U s ■ det Vi , Q' = U S Q Uj , and 
S' = U S SV? . U s and Vi are the SU(d) parts of 
the local unitary operations: U s = y/det U s U s and 
Vi = v^det Vi V . Therefore, we identify the transforma- 
tion in three sectors of the matrix structure: an explicit 
phase transformation <fi t— > <p' , a transformation closed in 
the space of positive definite Hermitian matrices Q ^ Q' , 
and a transformation S H> S' closed in SU(d). Since the 
roots of a complex number are multivalued functions, 
some caution is necessary in defining the above quanti- 
ties. For time varying unitary operations, we shall as- 
sume that any one of the possible roots is taken at the 
initial time, then the subsequent values must form a con- 
tinuous evolution as a function of time, so that <p(t) is a 
smooth, everywhere differentiable function. 



Following [18} , we shall define as cyclic those evolutions 
for which the final state is physically equivalent to the ini- 
tial one, in other words, the initial and final state vectors 
are related only by a global phase factor: a' = e l0 a . 
The geometric phase acquired by a time evolving quan- 
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turn state a(t) = e 1 ^^ 1 Q(t) S(t) is always defined as 



<t> g = arg(V>(0)|V(i))+t / dt (VWIV-W) (1) 
= arg{Tr[a t (0)a(t)]} +i J dt Tr [a + (t)d(f>] , 

that corresponds to the total phase minus the dynam- 
ical phase. At this point we can profit from the po- 
lar decomposition by investigating the contribution com- 
ing from each sector of the coefficient matrix. Let us 
consider a cyclic evolution over the time interval T : 
a(T) = e l6 a(0) . First, we identify the trivial phase evo- 
lution 4>(T) — 0(0) + A<fi . In the positive Hcrmitian sec- 
tor, if we write Q(T) — e teq Q(0), its defining condition 
imposes 8 q — , so that no phase contribution could stem 
from this sector. Finally, in the SU(d) sector we can have 



5(T) 



jo, 



5(0) , however, det 5(T) = e ld0 ° det 5(0), and 



since both 5(T) and 5(0) are SU{d) matrices, we arrive 
at 



2nn 
~d~ 



(2) 



where n = 0, l,2,...,d— 1. Therefore, only fractional 
phase values can arise from the SU(d) sector and the 
global phase acquired in a cyclic evolution is 



9 = A4>- 



2mr 



(3) 



Although SU(d) is simply connected, the topological na- 
ture of these fractional phases relies on the multiply con- 
nected manifold that represents the set of SU(d) matrices 
5 with the identification e l0E 5 = 5. 

We can also investigate the dynamical phase using the 
polar decomposition. At this point it will be useful to 
write: 



i<j>a(t) + e^ (t) \Q(t)S(t) + Q(t)S(t) 



(4) 



This expression can be used in the dynamical phase in- 
tegral appearing in eq.([2]). Using the normalization con- 
dition Tr[cr (t)a(t)] = 1, we see that the first term in 
eq. (0| will cancel out the trivial global phase Acj) , which 
do not contribute to the geometric phase, as would be 
expected. Then, by making o<t(t) = e"^*) 5* (t)Q(t) and 
using the cyclic property of the trace, one easily arrives 
at 
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+ i 



+ i 



dt Tr 



dt Tr 



Q(t)Q(t) 

(t)s(t)sHt) 



(5) 



The normalization condition implies in Tr[Q 2 (t)} = 1. 
Differentiating both sides with respect to t and using the 



cyclic property of the trace gives Tr Q(t)Q(t) 



0, so 



that the first integral in eq.® vanishes. Now, it is now 



important to identify the following invariants under local 
unitary evolutions: Tr[/jj ) ], p — l,...,d, where pj is the 
reduced density matrix of qudit j : 



p s = (a^a) T = (5tQ 2 5) 5 



and 



Q 2 



(6) 



(7) 



In fact, the invariants are j-independent since one eas- 
ily shows that Tr [pf] = Tr[pf ] = Tr[Q 2p ] . The first 
one (p = 1) is simply the norm of the state vector, as 
already stated. The second invariant is related to the 
I-concurrence of a two-qudit pure quantum state [l9j 

C = y/2(l — Tr[p 2 ]), so that its invariance expresses the 

well known fact that entanglement is not affected by local 
unitary operations. The p = d invariant can be rewritten 
in terms of the former and V = |deta|. In particular, 
for qubits we have C = 2 2?. For a given dimension d, 
the I-con currence ru ns between for product states and 
C m = \/2{d — l)/d for maximally entangled states. In 
order to exploit the role played by these invariants in 
the geometrical phase, we shall make them explicit in 
the expression of Q 2 in terms of the identity and the 
generators T n (n = 1,2,..., d 2 — 1) of SU(d) . As usual, 
the generators are traceless matrices normalized accord- 
ing to Tr [T n T m ] = 6 nm /2 . A general Q 2 matrix with 
Tr[Q 2 ] = 1 and Tr[Q 4 ] = 1 - C" 2 /2 can be written as 



Q\t) 



~d + ^ 



C 2 q(t) ■ T 



(8) 



where q is a unit vector in R . For maximally en- 
tangled states (C = C m ), the reduced density matrices 
become p s = pi = 1/d. 

The expression ([5]) can be used in the last term of 
eq. ([5]) to establish a connection between the geomet- 
rical phase and the concurrence of the two-qudit state 
undergoing a cyclic evolution. However, it will be useful 
to work out a more convenient parametrization of 5 5^ . 
First we recall that for a general invertible matrix A we 
have (H: 



d (det i 
~dT 



= det A Tr 



dt 



(9) 



Since the evolution of S(t) is closed in SU (d) (det S(t) 



1), we easily arrive at Tr 



5 



. Besides, note that 



S(t)S*(t) = -S(t)SHt), what allows us to define the 
velocity vector s € R d _1 such that 5 £t = is-T . Then, 
using the normalization condition for the generators and 
the fact that they are traceless matrices, we arrive at a 
compact expression for the geometrical phase: 



2nir 
~d~ 



C 2 <h q • dx , 



(10) 



where we defined dx — sdt . This key result shows that 
the geometrical phase for a cyclic evolution is composed 
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by a fractional contribution of topological nature and an 
integral contribution which depends on the history of 
the quantum state evolution, parametrized in the Her- 
mitian and SU(d) sectors of the polar decomposition. 
This integral contribution is weighted by entanglement 
and vanishes completely for maximally entangled states, 
for which only the fractional values are allowed. We next 
turn to a possible physical implementation of the frac- 
tional phases using entangled qudits encoded on spatial 
degrees of freedom of photon pairs generated by sponta- 
neous parametric down conversion (SPDC). 



III. POLARIZATION CONTROLLED 
TWO-PHOTON INTERFERENCE 

Let us consider the setup sketched in Fig.([T]). Using 
multiple slit masks, spatial qudits are encoded on the 
transversal paths of the photon pairs generated by SPDC. 
Beyond the slits, the wave field can be written as a super- 
position of the field transmitted through different slits. 
We shall refer to these contributions as the slit modes 
rj m (r) , where m is the slit index. The slit mode functions 
satisfy the orthonormality condition: 

/ d 2 r n* m {r) Vn {r) = 5 mn . (11) 

The positive frequency components of signal and idler 
vector field operators are 

E+ = E+ H e H + E+ V e v , 

E+ = E+ H e H + E+e v , (12) 

where e M (/i = H, V) is the horizontal (H) or vertical (V) 
polarization unit vector. Each polarization component is 
expanded in terms of the slit modes as 

v 

E+ = , (13) 

where greek indices refer to polarization modes and ro- 
man ones refer to slit modes. The annihilation operators 
a pfJl and b qv act on Fock states as usual: 

a pil b q „\ma, ne) \vac) , (14) 

where \vac) is the vacuum state and \ma,ne) is a two- 
photon Fock state corresponding to one signal photon 
passing through slit m with polarization a and one idler 
photon passing through slit n with polarization e. 

For type I phase matching and a vertically polarized 
pump beam, SPDC generates horizontally polarized pho- 
ton pairs. Meanwhile, the spatial correlations between 
signal and idler can be tailored by playing with the pump 
angular spectrum [l3T ] using suitable lenses (omitted in 




FIG. 1: Experimental proposal. NLC: Nonlinear crystal, A/2 
half waveplate, PBS: Polarizing beam splitter, SLM: Spatial 
light modulator, POL: polarizer, 8j (j — i,s) axe path phases 
added by the mirrors coupled to PZTs and D; (I = 1,2) are 
single photon detectors. 

the setup for simplicity) on the pump laser. These spa- 
tial correlations determine the two-photon path correla- 
tion through the slits. Therefore, the two-photon state 
after the slits can be written as 

IV'o) = ^ a m n\mH, nH) , (15) 

m,n 

where a m n is the probability amplitude of having one 
signal photon passing through slit m and one idler photon 
passing through n. Right after the slits, two half wave 
plates (A/2) rotate the photons polarization by 45°, and 
the two-photon quantum state becomes 

1^) = ( \mH, nH) + \mH, nV) 

+ \mV,nH) + \mV,nV)\ . (16) 

Then, each photon passes through a polarization con- 
trolled unitary gate composed by a Mach-Zehnder inter- 
ferometer with input and output polarizing beam split- 
ters (PBS), and a spatial phase modulator (SLM) in- 
serted on the vertical polarization arm. The signal and 
idler evolutions U s and V% are then implemented on 
the vertical polarization component of the spatial qudits 
encoded on signal and idler photons, respectively. Of 
course, the efficient implementation of the unitary oper- 
ations require that the slits are imaged on the SLM's. For 
the sake of simplicity, we shall omit the lenses required 
for this imaging scheme. After the controlled gates, the 
two-photon quantum state becomes 

|^ 2 ) = ^ ( e <( " s+9i) \mH, nH) + e w °V t \mH, nV) 

+ e t6 'U s \mV, nH) + U s V t \mV, nV)} , (17) 

where S and 9i are longitudinal phases added to signal 
and idler when they cross the Mach-Zehnder interferom- 
eters. The displacement of a piezoelectric ceramic (PZT) 
coupled to one of the Mach-Zenhder mirrors can do this 
task. 
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Finally, signal and idler are sent to different input ports 
of the same PBS so that only the components \nH, raH) 
and \nV, mV) will contribute to coincidences. Two po- 
larizers (POL) oriented at 45° are placed before detectors 
in order to erase the polarization information. Then, the 
detected field operators are: 



Et = -±={E+ H + iE+) , 

and the coincidence count is proportional to: 

(ip 2 \E^EzE+E+\i> 2 ) = \\E+E+\^ 2 )\\ 2 



(18) 



(19) 



where E A = (E^Y and j = 1,2. From eq. ([T4|) . one eas- 
ily sees that the \mH,nV) and \mV,nH) contributions 
vanish and the normalized coincidence function reduces 
to 



C(ri , r 2 ) 



e W a™ E+ H E+ H \mH, nH) 



a' mn E+ v E+ v \mV,nV) 



(20) 



where 8 = 8 S + 6% — tt . The coefficients a' mn result from 
the local unitary transformations on the two-qudit state: 



mV, nV) 



\mV,nV) . (21) 



The two-photon interference pattern is 
1 



C(ri , r 2 ) 



^77 m (ri)77„(r 2 ) (a r , 



(22) 



We are interested only in the phase acquired by the 
qudits when they evolve under the action of the uni- 
tary gates U s and VI. Since the slit modes are orthonor- 
mal, the integrated coincidence count eliminates the spa- 
tial interference between different slits, so that only the 
Mach-Zehnder interference shows up in the two-photon 
correlations. This corresponds to large aperture detec- 
tion, insensitive to the detailed spatial structure of the 
two-photon quantum correlations. However, the Hilbert 
space dimension remains manifested through the two- 
qudit coefficients in the integrated coincidence function: 



C = J d 2 v± d 2 r 2 C(n , r 2 ) 

= It. 



(23) 



Note that polarization has been used as a subsidiary 
degree of freedom to provide two paths for the evolution 



of the spatial qudits. In fact, the interference described 
by ea. (f2"3"f can be thought of as resulting from the su- 
perposition between the initial two-qudit state \ipo) = 
J2mn a mn| mn ) an d the evolved state \<p) = U a Vi\tpo) such 
that 



C = 



\<Po) + \<P) 



V2 

| {<p \tp)\ cos 



(0-7) 



(24) 



where 7 = axg((po\<p) . The visibility of the interference 
fringes is the absolute value of the state overlap (ifio\f) , 
while the interference phase is the overlap argument. For 
a cyclic evolution, \ip) — e n \tpo), the interference recovers 
maximal visibility with the fringes shifted by 7 . Finally, 
the topological phases can be implemented by continu- 
ously varying the unitary gates according to local time 
dependent parameters, adjusted by the signal and idler 
SLM's respectively. Next, we give some numerical ex- 
amples of the possible operations and the expected in- 
terference patterns that can reveal both the fractional 
topological phases and the role of entanglement. 



IV. NUMERICAL RESULTS 

We can illustrate the fractional phases with the sim- 
plest example of two qutrits (d = 3) under the action 
of diagonal unitary operations U s 



and Vi = diag[e 



diag[e^\e^ 2 ,e l<te ] 
e* X3 ] , where <fi n and \n are phase 
factors imposed by the spatial light modulators (SLM) 
to the photon path states. In order to guarantee that 
no trivial dynamical phase is added by the operations, 
we assume that J2 n = J2n Xn — , in other words, 
we restrict the two-qutrit evolutions to SU(3) ® SU(3) 
operations. 

The coefficient matrix evolves as a' — U s a , so that 



a,. 



(25) 



and the coincidence count becomes 



C 



zE ^ 



e ie + e i[0 m+x „] 



E 



a mn cos 



<t>rn + Xr. 



(26) 



In principle, the two SLMs can be operated indepen- 
dently and implement totally different local evolutions. 
However, for simplicity, we shall first consider an iden- 
tical operation capable to evidence both the fractional 
phases and the signatures of entanglement. Our cyclic 
operations will be decribed by a real parameter t in the 
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(a) 
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FIG. 2: Coincidence interference for (a) a maximally entan- 
gled state and (b) a product state, corresponding to the uni- 
tary operations given by eqs. (l27)) . for three values of t: t = 
(continuous line), f = 0.5 (dashed line) and t = 1 (blue line). 



interval t e [0, 1] . As a first example, let us consider 



Xi 



2t-(2t-l) H [t 



2tt 

92 = X2 = - Y f ' 

^3 = X3 = | (2i - 1) H (t - \ 



(27) 



where H(t) is the Heaviside function. Eqs. (|27l) decribe 
a continuous phase evolution in each qudit component, 
corresponding to a local SU(d = 3) operation indepen- 
dently applied to each qudit. Under this SU(3) <X> SU(3) 
operation, we may compare a maximally entangled state: 



We 



1 

71 



(|00) + |11> + |22» , 



(28) 



(ctmn = fimn/VS) with a product state having the same 
single qudit population distribution: 



I^P> = | (|0) 



ID 



|2))®(|0) + |1) + |2)), (29) 



(a mn = 1/3). The resulting coincidence interference 
given by eq. (l2l)l) is shown in Figs. [5^ and ^jp. As we 



can see, the behaviour is considerably different for prod- 
uct and entangled states. For the entangled state, the 
interference fringes completely disappear when t = 1/2, 
and the two-qudit entangled state becomes orthogonal 
to the initial state. Then, the fringes reappear with the 
fractional phase shift 2-7r/3, reaching maximal visibility 
again at t = 1. On the other hand, the product state ex- 
hibits a remarkable different behaviour, the fringes never 
disappear, and maximal visibility is not recovered when 
t = l. 

Let us now consider a different kind of evolution con- 
trolled by independent parameters t s and tf. 

7T 

(Pi = 02 = ^ IS , 

2tt 



Xi = X2 = g U , 
2tt 

X3 = -y*! ■ 



In this case, the coincidence is 



C e = - cos 2 



1 2 
-COS 



2£tt 



(30) 



(31) 



for the entangled state of eq. (f28f , and 



tn 

y 



2tTT 



1 + cos (ttt) cos 



tTT 



(32) 



for the product state of eq. ((29|) . where we defined t = 
(t s + ti)/2 and r = (t s — ti)/2 . As we can see, the coin- 
cidence interference is completely different in each case. 
It does not depend on the relative parameter t for the 
entangled input, while it does for the product one. The 
behaviour of the fringe visibility is also considerably dif- 
ferent, even for t s = t. L = t (r = 0) . In Fig. [3] the 
coincidence interference is plotted for r — and different 
values of t . In this case, the evolved entangled state never 
becomes orthogonal to the initial one, but the fringes re- 
cover maximum visibility for t = 1, with the fractional 
phase shift 27r/3. However, the same operation does not 
take the product state through a cyclic evolution and the 
associated interference never recovers maximum visibil- 
ity. 

Finally, in order to observe the effect of the Hilbert 
space dimension, it is instructive to consider a situation 
in which only two components of the qutrits are operated, 
and compare their evolution with a genuine pair of qubits. 
For example, let us make: 



9i 

03 

Xi 

X3 



1 n 

h = 7T t s , 



= 



7T 

-X2 = U , 



(33) 
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(a) 
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single qubit population distribution: 



1 



|^) = -(|0) + |1))®(|0) + |1» - 



(«mn = 1/2): evolving under 



91 



TT 

-02 = ^ t S 
TT 



In this case, the coincidence is 



(37) 



(38) 



(b) 
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FIG. 3: Coincidence interference for (a) a maximally entan- 
gled state and (b) a product state, corresponding to the uni- 
tary operations given by eqs. (f30)) . for three values of t: t = 
(continuous line), t = 0.5 (dashed line) and t = 1 (blue line). 



In this case, the coincidence is 



2 



cos ( 



6 



for the entangled state 



[1 + 2cos(?ri)] 



and 



(34) 



Co 



cos 9 



+ COS (-7T t) + COS (tT t) 



2 cos (| t) cos (| 



(35) 



for the product state (|29j) . As before, the coincidence 
count associated to an entangled input is independent 
of the relative control parameter r, while the same does 
not hold true for the product input. However, as can be 
readily seen in Figs. 0Ji and0]3, neither the entangled nor 
the product inputs are driven through a cyclic evolution 
in the full range of the total control parameter t . The 
interference fringes never recover maximal visibility. 

We can compare these results with those expected for 
a pair of qubits (d = 2) prepared either in the entangled 
state 



1 



We 



V=2 m 



|H)) 



(36) 



(a r , 



fimn/ \/2) or in the product state having the same 



C e (d = 2) = - [1 + cos (tt t) cos 9] , 



(39) 



for the entangled state (|36|) . and 



C p (d- 



cos 9 



cos (tt t) + cos (ttt)] , (40) 



for the product state (|3"T|) . These coincidence counts are 
shown in Figs. [5^ and[5}) for r = and the full range oit . 
Now, the entangled state is taken through a cyclic evolu- 
tion and recovers maximal visibility with the topological 
phase tt, as expected for maximally entangled qubits 0- 
llOj . On the other hand, the product state does not com- 
plete a cyclic evolution, it ends up in a state orthogonal 
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FIG. 4: Coincidence interference for (a) a maximally entan- 
gled state and (b) a product state, corresponding to the uni- 
tary operations given by eqs. (|33[) . for three values of t: t = 
(continuous line), t = 0.5 (dashed line) and t = 1 (blue line). 
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sion of the qudits being operated. This leads to a natural 
conjecture that the quantities being measured are related 
to entanglement and dimensionality witnesses, although 
we shall leave the formal investigation of this aspect to a 
future work. 



V. CONCLUSIONS 

As is well-known, qudit gates based on topological 
phases are a potentially robust means for implementing 
quantum algorithms [2l| - |23j . In this work, we proposed a 
physical implementation of fractional topological phases 
acquired by entangled qudits operated by local unitary 
transformations. Spatial qudits encoded on quantum cor- 
related photon pairs are a suitable framework for this 
investigation, since the qudit components may be effi- 
ciently addressed by spatial light modulators. The ex- 
pected coincidence interferences can reveal both, the role 
played by entanglement as well as the Hilbert space di- 
mension. Then, a natural conjecture arises as to whether 
the measured quantities are related to entanglement and 
dimensionality witnesses, what motivates further investi- 
gation. 
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FIG. 5: Coincidence interference for a (a) maximally entan- 
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